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DISCOVERY OF THE MOST ISOLATED GLOBULAR CLUSTER IN THE LOCAL UNIVERSE 
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ABSTRACT 

We report the discovery of two new globular clusters in the remote halos of M81 and M82 in the 
M81 Group based on the Hubble Space Telescope archive images. They are brighter than typical 
globular clusters (My = —9.34 mag for GC-1 and My = —10.51 mag for GC-2), and much larger 
than known globular clusters with similar luminosity in the Milky Way Galaxy and M81. Radial 
surface brightness profiles for GC-1 and GC-2 do not show any feature for tidal truncation in the 
outer part. They are located much farther from either of M81 and M82 in the sky, compared with 
previously known star clusters in these galaxies. Color-magnitude diagrams of resolved stars in each 
cluster show a well-defined red giant branch (RGB), indicating that they are metal-poor and old. We 
derive a low metallicity with [Fe/H] rj —2.3 and an old age ~ 14 Gyr for GC-2 from the analysis of the 
absorption lines in its spectrum in the Sloan Digital Sky Survey in comparison with the simple stellar 
population models. The /-band magnitude of the tip of the RGB for GC-2 is 0.26 mag fainter than 
that for the halo stars in the same field, showing that GC-2 is ^400 kpc behind the M81 halo along 
our line of sight. The deprojected distance to GC-2 from M81 is much larger than any other known 
globular clusters in the local universe. This shows that GC-2 is the most isolated globular cluster in 
the local universe. 

Subject headings: galaxies: evolution — galaxies: individual (M81, M82) — galaxies: star clusters: 
general — galaxies: groups: individual (M81 Group) 



1. INTRODUCTION 

According to the current paradigm of large scale struc- 
ture formation, individual galaxies, galaxy clusters, and 
groups are formed via hierarchical merging of galaxies. 
Globular clusters are a powerful tool to test this hypoth- 
esis. Globular clusters are often found in and around 
galaxies. With the advent of wide field surveys, globu- 
lar clusters are sometimes found in the remote halo of 
galaxies in the Local Group. A small number of globu- 
lar clusters are found beyond 30 kpc from the center of 
the Milky Way Galaxy (MWG) and M31 (jHarrisI 119961 : 
iHuxor et al"1l2008t iGalleti et aTIl2007ft . while most of the 
globular clusters are much closer to the galaxy center. 
A few clusters are found also in the remote halo of less 
massive ga laxies such as NGC 6 8 22 (iHwang et al.|l20ilft 
and M33 (jStonkute et all l2008t iCockcroft et all 12011ft . 
Recently intracluster globular cl usters are also found 
in nearby galax y clusters : Virg o (Willia ms et al.l [20071 : 
iLee et all 12010ft . Coma (jPeng et al.l 12011ft . and Abell 
1835 (|West et al.l 12011ft . Globular clusters in galaxies 
provide fossil records for early formation of spheroidal 
components in the collapsing phase, while globular clus- 
ters far from galaxies reveal clues for later growth of 
galaxies via accretion. 

The M81 Group, one of the nearest galaxy groups, is an 
excellent laboratory for studying the property of dwarf 
galaxies and star clusters as well as starburst g alaxies 
and intergalactic medium (jChiboucas et al.l 120 09!) . The 
main galaxy located in the center of t he group is M81, 
surrounded by 29 member galaxies (jChiboucas et al.l 
20091: iMakarov fe~K arachcntscv 12011ft . According to 
Makarov fc KarachentsevT^Oll ). the systemic velocity 



with respect to the Local Group and the velocity dis- 



persion of the M81 Group are vlg = 193 km s _1 and 
a v = 138 km s _1 , respectively, and the total mass derived 
from the velocity of the members is M — 3.89 x 10 12 M Q . 
The M81 Group includes also a famous s tarburst galaxy 
M82. iKarachentsev fc K ashibadzc (2006) estimate from 
kinematics and distances of the member galaxies that 
M81 is twice as massive as M82 and the mass of the M81 
Group is 77% of the mass of the Local Group. 

The distance to M81 is known to be 3.63 ± 0.14 Mpc 
((m-M)o = 27.80±0.08) derived using the tip of the red 
giant branch (TRG B) method from deep HST images as 
well as Cepheids bv iDurrell et al.l ([2010). The distance 
to M81 is estimated to be similar to that to M82 based on 
the same TR GB method, 3.55 ± 0.11 Mpc ((to - M) 
27.7'i ■ 0.07 : (|Lee fc Lim 112012ft . showing that these two 
galaxies are at a similar distance from us. 

Previous studies found numerous star clusters in M81 
and M82, whic h are mostly located in the main body 
of each galaxy flMayva et al.l 120081: iNantais et al.l 1 2010 1: 
i Nantais fc Huchral l2010t iSantiago-Cortes et al.l 120 lOL 
120111 : INantais et al.l 120111 ). We have been searching for 
globular clusters in a remote halo region of each galaxy, 
much farther from previous surveys. In this paper we 
present a discovery of two globular clusters in the re- 
mote region of M81 and M82. This paper is composed as 
follows. In Section 2, we describe data used and the glob- 
ular cluster search method. Section 3 presents the dis- 
covery of two new globular clusters, the color-magnitude 
diagrams of resolved stars in each cluster, and distance 
measurements. We also present the measurements of 
metallicity, age, and [a/Fe] using the spectrum of one 
globular cluster. Finally we derive the surface brightness 
profiles of the new globular clusters. Primary results are 
discussed and summarized in the final section. 
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Fig. 1. — (a) A gray scale map of the digitized sky survey for the 
central region of the M81 Group, showing the positions of GC-1 
and GC-2 discovered in this study. North is up and east to the 
left, (b) and (c) Gray scale maps of the F81AW images for GC-1 
and GC-2. 

2. DATA AND CLUSTER SEARCH 

We used the images taken with the Hubble Space Tele- 
scope (HST)/ Advanced Camera for Surveys (ACS) and 
Wide Field Camera 3 (WFC3) in the archive to search 
for globular clusters in the remote halo of M81 and M82. 
There are ~50 sets of multi-band images for the remote 
regions in the M81 group available. The images in which 
we found new globular clusters are ACS/WFC F606W 
and F814W images (obtained with exposure times 850 s 
and 690 s, respectively), and WFC3/UVIS F606W and 
F814W images (obtained with exposure times 735 s and 
1225 s, respectively). Both fields come from the HST 
program 11613 (P.I. : de Jong). 

We expected that globular clusters far from galaxies 
are more extended than those close to galaxy centers, so 
that some stars in the outer region of these globular clus- 
ters at the distance of the M81 Group can be resolved 
in the HST images. We searched all these fields visually 
with a primary aim to find globular clusters. To identify 
an object as a globular cluster, we used the following vi- 
sual criteria : (1) round appearance, (2) lack of a smooth 
extended halo typical of background elliptical galaxies, 
and (3) being surrounded by an excess of resolved stars. 

3. RESULTS 

3.1. Discovery of New Globular Clusters 

Through visual search we discovered two new glob- 
ular clusters, JM81GC-1 and JM81GC-2 (called GC-1 
and GC-2 hereafter), in ACS/WFC and WFC3/UVIS 
fields, respectively. Figure [Ha) shows the location of 
these globular clusters. GC-1 is 15'. 90 west of M82 and 
29'. 58 north of M81 (corresponding projected distances 
are 17.07 kpc and 31.77 kpc, respectively). GC-2 is 13'. 34 
west of M82 and 37'. 25 north of M81 (corresponding 
projected distances are 14.32 kpc and 40.00 kpc, respec- 
tively). Their positions are RA(2000)=09 ?l 53™ 26.22 s , 



Dec(2000)=69° 31' 17". 5 for GC-1, and RA(2000)=09 ?l 
53 m 20.17 s , Dec(2000)=69° 39' 16".4 for GC-2. 

We checked the images o f these clusters i n the Sloan 
Digital Sky Survey (SDSS) (lYork et al.ll2000l) . GC-1 and 
GC-2 appear as extended source in the SDSS images and 
were classified as galaxies (SDSS ID : J095326+693117 
for GC-1 and J095320+693916 for GC-2). However, 
F814W images in Figure 1(b) and (c) show some re- 
solved stars in the outer region of these clusters, proving 
that they are genuine star clusters. We checked the ex- 
istence of any nearby d warf galaxies around th ese two 
clusters using the data in lChiboucas et ail (|2009D as well 
as the HST images, but found none. 

3.2. Color-Magnitude Diagrams of Resolved Stars 

We derived instrumental magnitudes of point sources 
in the images using the IRAF/ DAOPHOT package that 
is design ed for point sp read function (PSF) fitting pho- 
tometry ([Stetson! [T994) . We used 2-cr as the detection 
threshold, and derived the PSFs using isolated bright 
stars in the images. We applied aperture correction 
derived from isolated bright stars. We calibrated 
the instrumental magnitudes to the Vega magnitude 
syste m using photome tric zeropoints for A CS /WFC 
( http: / /www. stsci.edu/hst /acs /analysis /zeropoints / # 
t ablest art ) and W FC3/ U V IS ( http : / / www . stsc i.edu/] 
hst/wfc3/phot_zp_lbn). Then we convert ed this system 
to Jo hnson-Cousins VI system using iSirianni et al.l 
<f2005h - 

Figures Ufa)— (d) display the color- magnitude dia- 
grams for GC-1 and GC-2 as well as corresponding fields 
(called Field 1 and Field 2, respectively). Because the 
crowding is severe in the central region of the globular 
clusters, we plotted the stars at 0".8 < r < 6".0 from 
the cluster center. The aperture with radius of 6" covers 
about 96 % of the total luminosity of each star cluster. In 
the case of fields, we plotted the stars at 15" < r < 100". 
Both clusters show a relatively well-defined red giant 
branch (RGB), indicating that they may be old globular 
clusters. 

3.3. Distance Estimation 

We derived the distance to these clusters as well 
as to the corre sponding halos using the TRGB 
meth o d (ILee et al.lll993tlSakai et al.lll996t lMendez et all 
| 2002t IMcConnachie et all 120041: IMouhcine et all 120101: 
iConn et al.l 1201 lTT In Figure e) the /-band luminos- 
ity function of the red giants in GC-1 and Field- 1 shows 
a sudden jump at J « 24.0, which corresponds to the 
TRGB. However, the TRGB magnitude for GC-2 is 
about 0.3 mag fainter than that for Field-2 in Figure 
HJg), showing that GC-2 may be behind the halo stars. 

Using the edge-detecting algorithm, we determined the 
TRGB magnitude more quantitatively. We calculated an 
edge-detection response function E{m) (— $(m + <r m ) — 
$(m— o~ m ) where $(m) is the luminosity function of mag- 
nitude m and a m is the mean photometric error within a 
bin of ±0.05 mag about magnitude m), and we weighted 
it according to the Poisson noise of the luminosity func- 
tion E(m) y '<&(m) (Mc ndez et al.ll2002T) . as shown in Fig- 
ures Hp) and (h). 

Thus derived TRGB magnitudes are 7 T rgb = 24. 039 ± 
0.021 for GC-1, 24.069 ±0.011 for Field- 1, 24.295 ±0.036 
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Fig. 2. — (a) — (d) I-(V — I) color-magnitude diagrams of the stars in the region at 0".8 < r < 6".0 for GC-1 and GC-2 and in their 
corresponding fields at 15" < r < 100". The curved lines represent the loci of the RGB in the Milky Way globular clusters with a range 
of metallicity ([Fe/H] = -2.17, -1.58, —1.29, and —0.71, respectively), shifted according to the foreground reddening and derived distance. 
Errorbars represent the mean errors, (e)-(h) /-band luminosity functions of red giants with 0.7 < (V — I) < 1.9 (e and g), and weighted 
edge-detection responses (f and h). Arrows represent the positions of the TRGB. 



for GC-2, and 24.036 ± 0.022 for Field-2. The errors 
for the TRGB magnitudes were determined using boot- 
strap resampling method with one million simulations. 
In each simulation, we resampled randomly the RGB 
sample with replacement to make a new sample of the 
same size. We estimated the TRGB magnitude for each 
simulation using the same procedure and derived the 
standard deviation of the estimated TRGB magnitudes. 

The mean color of the TRGB is derived from the colors 
of the bright red giants close to the TRGB and is cor- 
rected for the foreground reddening (E( B — V) = 0.09 
for G C-1 and E{B - V) = 0.10 for GC-2 (jSchlegel et al.l 
1998)): (V — I)o trgb = 1.32±0.02for GC-1, 1.41±0.03 
for Field-1, 1.31 ± 0.04 for GC-2, and 1.49 ± 0.04 for 
Field-2. We derive the intrinsic J-band magnitude of the 



TRGB using the relation between the bolometric mag- 
nitude Afbol and the bolometric correction BC : M\ 
= Mboi — BC\ . We calculate the bolometric magni- 
tude using Mboi = — 0.19[Fe/H] — 3.81 and the bolomet- 
ric correction using BQ n = 0.8 81 — 0.243(V — I)o trgb 
(|Da Costa fc Armandron1[T990l ). [Fe/H] is derived from 
the mean color of the RGB stars. 

The mean color of the RGB stars 0.5 mag fainter than 
the TRGB is derived from the colors of the bright red 
giants close to this magnitude: (V— J)o -35 = 1.21±0.02 
for GC-1, 1.22 ± 0.03 for Field-1, 1.21 ± 0.03 for GC-2, 
and 1.31 ±0.02 for Field-2. From these we derive [Fe/H]: 
[Fe/H] = -2.23 ±0.11 for GC-1, -2.18 ±0.13 for Field- 
1, -2.23 ± 0.12 for GC-2, and -1.84 ± 0.10 for Field-2. 
Using these we obtain M^trgb = —3.95 for GC-1, —3.93 
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Fig. 3. — (a) Lick line indices for H/3 versus [MgFe]' for GC-2, useful for measuring the me tallicity and age (dots with errorbars). Grids 
represent the SSP models for various values of [Z/H] and age[Gyr] at [a/Fe] = -0.05 given by Thomas ct al. (2003). (b) <Fe> versus Mg2 
useful for measuring [a/Fe]. The grids are for an age 14.9 Gyr. 



for Field-1, -3.95 for GC-2, and -3.98 for Field- 2. 
Thus we finally calculate the distance modulus using 
(m — M)o = /o trgb — Mi trgb ■ The distances de- 
rived are (m - M) = 27.80 ± 0.03 (d = 3.63 ± 0.05 
Mpc) for GC-1, 27.81 ± 0.03 (d = 3.65 ± 0.05 Mpc) for 
Field-1, 28.04±0.04 (d = 4.05±0.08 Mpc) for GC-2, and 
27.81±0.03 (d = 3.65±0.05 Mpc) for Field-2. This value 
is sim ilar to the previous estimates for M81, 3.63 ± 0.14 
Mpc (IDurrell et ail \2(M IGerke et al.l IMtl and M82, 
3.55 ± 0.11 Mpc (|Lee fc Lim 1 12012ft. While GC-1 is at 
the same distance as M81, GC-2 is about 400 kpc behind 
the M81 halo along our line of sight. 

3.4. Spectral Line Analysis for GC-2 

GC-2 was previously observed and classified as a 
galaxy in the SDSS and its optical spectrum is avail- 
able in the SDSS. The spectrum of GC-2 shows several 
absorption lines typical for globular clusters. We derived 
[Fc/H], [a/Fe], and age for GC-2 from the comparison of 
Lick line index diagra ms with the simple s tellar popula- 
tion (SSP) models bv lThomas eTaTl (|2003ft . 

Figure [3] displays H/3 versus [MgFe]' diagram, and 
<Fe> versus Mg2 diagram for GC-2. The com- 
posite index [MgFe]' is a combination of magne- 
sium and iron-sensitive indices, defined as [MgFe]' = 

v/Mgb(0.72 • Fe5270 + 0.28 • Fe5335). It is an excellent 
metallicity indicator, because it is completely indepen- 
dent of [a/Fe], and this behavio r is almost independ ent of 
the adopted age or metallicity (jThomas et~aT1l200l . H/3 
is an efficient indicator for age. We used the line index 
data for these clusters provided by the SPSS . The values 
derive d foll owing the technique d e scribe d in lPuzia et al.l 
(2005) and lPark. Lee, fc Hwang I d20H are [Fe/H] = - 
2.3 ± 0.12, [a/Fe] = -0.05 ± 0.40, and age = 14.9 ± 1.0 
Gyr. This metallicity is consistent with the value de- 
rived from the color of the RGB. These show that GC-2 
is indeed very metal-poor and old. 

3.5. Spectral Energy Distribution Fit 

We also derived age and mass from the spectral energy 
distribution fit using SDSS ugriz magnitudes in com- 



parison with the SSP model given by lBruzual fc Charlotl 
([2003ft (u = 20.50 ± 0.12, g = 19.13 ± 0.02, r = 
18.50±0.02, i = 18.21±0.02, and z = 18.03±0.05 for GC- 
1, and u = 19.51±0.06, g = 18.25±0.01, r = 17.60±0.01, 
i = 17.31 ± 0.01, and z = 17.10 ± 0.03 for GC-2). De- 
rived ages are log(age[y]) ~ 10.2 for both clusters as long 
as we adopt Z = 0.0001, and derived cluster masses are 
log(M/M Q ) ~ 6.40 for GC-1, and 6.85 for GC-2. This 
age for GC-2 is also consistent with the value derived 
from the spectrum analysis. Thus these clusters are very 
massive and old. 

3.6. Surface Photometry 

We derived surface photometry of GC-1 and GC-2 us- 
ing IRAF/ELLIPSE task from the HST images. Figure 
4 displays the radial profiles of the surface brightness for 
F81AW images. The radial profiles of the inner region 
look similar to the King profiles. However, the outer 
parts in the radial profiles do not show any tidal cut- 
off, but continue to decrease smoothly. From the King 
model fits for the inner region at 0".0 < r < 4".0, we 
derived core radii (r c ), 0".0755 ± 0".0003 (1.329 ± 0.005 
pc) for GC-1 and 0".1269 ± 0".0003 (2.491 ± 0.006 pc) 
for GC-2. Also we fit the data usin g the model with a 
power law form in the outer region in lElson et al.1 (|1987ft . 

,,2 -7/2 

XX r p) — So (1 + ^) where the scale radius a is re- 
lated with the King core radius as r c — a\j2 2 h — 1 for 
r t > r c . There is a break at r = 2".80 (55 pc) for GC-2 
so that we fit the data in two parts: a = 1.96, 7 = 2.48 
for GC-1, and a = 3.16, 7 = 2.42, a = 3.27, 7 = 3.62 for 
GC-2. The radial surface brightness profiles of GC-1 and 
GC-2 gradually decrease out to r = 15". By integrating 
the radial surface brightness profiles to r = 15", we de- 
rived integrated magnitudes of the star clusters, which 
we shall present in Section 4.1. From these we derived 
half light radii (r h ) : 0".352 ± 0".006 (6.13 ±0.01 pc) for 
GC-1 and 0".511 ± 0".005 (9.81 ± 0.01 pc) for GC-2. 

4. DISCUSSION AND SUMMARY 
4.1. Luminosity- Size Relation 
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Fig. 4. — F814Vy-band surface brightness profiles for GC-1 (a) and GC-2 (b) (dots with errorbars). The thick solid lines represent the 
King model fit and the dotted and dashed lines represent the EFF power-law model fit. 
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Fig. 5. — Half light radii (r h [pc]) versus My for GC-1 
and GC-2 (star symbols) in comparison wi th globular c lusters 
in the Milky Way Galax y (open circles') jHarrisI 119961) . M31 
(filled cir cles and crosses'! IIHuxor et al.l 120091 . l20TTf t. M81 (open 
squares) UNantais et al.l 120101 . 1201 ll) , and M82 (open triangles) 
IILim. Hwang, fc Lee 1120121) . 



Integrated magnitudes and colors of GC-1 and GC- 
2 are derived, respectively, Vq — 18.46 ± 0.01 and 
(V - I) = 0.91 ± 0.01 and V = 17.53 ± 0.01 and 
(V — 7)o = 0.90 ± 0.01. Corresponding absolute magni- 
tudes are M v = -9.34±0.01 for GC-1 and -10.51±0.01 
for GC-2, derived using our TRGB distances to the 
clusters. In Figure [5] we plot half light radii (r^) ver- 
sus absolute magnitude My for GC-1 and GC-2 in 
comparis on with glob ular cl usters in the Milky Way 
Galax y (lHarrisI [l996h. M81 dNantais et al.l [201lh and 
M82 (|Lim. Hwang, fc Lee I [2012D . We also plotted the 



data for globular clusters and extended star clusters in 
M31 dHuxor et al.l [20091 l20ll. 

Ivan d en Bcrgh & Mackcv (2004) presented a boundary 
relation between normal globular clusters and extended 
globular clusters: \og(r h ) = 0.2M V + 2.6. GC-1 and GC- 
2 are found to be located above this boundary line, like 
the case of u> Cen and NGC 2419 in our Galaxy, and May- 
all II in M31. Thus GC-1 and GC-2 are larger than typ- 
ical globular clusters with similar luminosity, and they 
are much brighter than the extended star clusters in M31. 
The latter globular clusters (w Cen, NGC 2419, and May- 
all II) show some peculiar features so that they are often 
considered to be remnants of dwarf galaxies. Therefore 
GC-1 and GC-2 also may be in the same vein. 

4.2. Intragroup Globular Clusters 

At the moment the most isolated globular cluster in 
the Local Group is known to be MGC1, located at the 
projected distance from M31, 117 kpc. Considering it 
is ~ 160 kpc closer than M31, its deproj ected distance 
from M31 was derived to be 200 ± 20 kpc (jMackev et al.l 
120101) . It is brighter (M v = -9.2) than typical globular 
clusters, but fainter than GC-1 in M81. In the Milky 
Way Galaxy, the most distant globular clusters and their 
galactocentric distances are NGC 2419 (91.5 kpc), Pal 3 
(95.9 kpc), Eridan us (95.2 kpc) , Pal 4 (111.8 kpc), and 
AM-1 (123.2 kpc) (IHarrislll996l) . 

On the other hand the most distant satellite dwarf 
galaxies in the two main Local Group galaxies are Leo I in 
the Milk y Way Galaxy and And XXVIII at d > 350 kpc 
in M31 ([Slater et al.l 1201 lT ) . Leo I is a dwarf spheroidal 
galaxy located at 270 kpc, being considered long as 
the most distan t satellite of the Milky Way Galaxy 
(|Lee et al.l fl993h . And XXVIII is a newly discovered 
dwarf galaxy about 100 kpc closer than M31. Its de- 
projected distance from M31 is estimated to be 365^ 7 
kpc (jSlater et al.ll201ll ). Thus the most remote globular 
clusters in our Galaxy and M31 are closer than the most 
distant satellite dwarf galaxies. This raises interesting 
questions: Can there be any globular clusters more dis- 
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tant than the most distant satellite galaxies in a galaxy? 
If so, what are they? 

GC-2 lies 406 ± 97 kpc behind M81 along our line of 
sight, according to our TRGB distance estimates to GC-2 
and Field-2. Considering this and the projected distance 
in the sky we derive a three-dimensional distance from 
M81, 408±97 kpc. Although it is closer to M82 than M81 
in the sky, its three dimension distance is significantly 
larger than the projected distance from M82 (also from 
M81). Also M81 is twice as massive as M82 and is the 
most massive member in the M81 Group. These indicate 
that the main host for GC-2 may be, if any, M81 rather 
than M82. Therefore GC-2 is the most isolated globular 
cluster among the known globular clusters in the local 
universe. In addition, it is more distant than any known 
satellite dwarf galaxies around M81. The radial velocity 
for GC-2 given in the SDSS is 159 ± 4 km s" 1 , which is 
much larger than the y alue for M81 (—35 ± 4 km s^ 1 ) 
(jChynoweth et all [2008) . Thus GC-2 is receding with a 
velocity w 200 km s _1 from M81. If this were moving at 
this velocity along the radial orbit, it would have taken 
about 2 Gyr to reach the current position from the center 
of M81. Thus GC-2 can be considered as an intragroup 



globular cluster wandering in the M81 Group. 

How can GC-2 be that far from M81, receding with 
a velocity of « 200 km s^ 1 ? The origin of GC-2 is not 
clear. Possible scenarios are as follows. First, it might 
have ejected during the interaction of three main galaxies 
about 2 Gyr ago. However, there are no evidence sup- 
porting this at the moment. Second, it may be one of 
the primordial globular clusters that formed early in iso- 
lation. If so, where are others? Third, it may a remnant 
of a dwarf galaxy that accreted to M81 and is receding 
now. If so, how could it survive during the perigalac- 
tic passage around M81? These possibilities need to be 
investigated with observations or simulations. 
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